
X-642-68-51 

PREPRINT 

NASA m h b?> II J 



) 


STELLAR OPACITY PROGRAM MANUAL 

"STOP” 

■ ■ PRELIMINARY' REPORT 


GPO PRICE $ 


CFSTI PRICE(S) $ 

A A THOMAS E. MICHELS 

Hard copy (HCI 0,0 ' IRAKASKEL 

Microfiche (MF) SHEILA MEYER 


ff 653 July 65 



FEBRUARY 1968 




GODDARD SPACE FLIGHT CENTER 

GREENBELT, MARYLAND 



STELLAR OPACITY PROGRAM MANUAL 


"STOP" 


Thomas E. Michels 
Ira Kaskel 
Sheila Meyer 


February 1968 


GODDARD SPACE FLIGHT CENTER 
Greenbeit, Maryland 


Pi<tCtuu\ 


cp,i^G PAC 


'E BLANK. NOT: FILMED. 


STELLAR OPACITY PROGRAM MANUAL 
"STOP" 


Thomas E. Michels 
Goddard Space Flight Center 
Greenbelt, Maryland 

Ira Kaskel 
Sheila Meyer 

Computer Applications, Inc. 
New York, New York 


ABSTRACT 

The new Stellar Opacity program, "STOP" is described. "STOP" is the 
FORTRAN version of the IBM 7090 FAP, "OPACITY" and "IONIC" programs 
developed by Dr. Arthur N. Cox at the Los Alamos Scientific Laboratory. A 
Rosseland Mean Continuous Opacity is computed and the following processes 
are included: Bound-free and free-free absorptions including those of the H - 
ion, Compton electron scattering, H 2 and H+H molecular absorptions, and Elec- 
tron conduction. Equations used for these processes as well as the "Mayer 
Independent Electron" and "Ionic" methods of computing occupation numbers 
and other equation of state data are discussed. "STOP" is available for general 
usage. 
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STELLAR OPACITY PROGRAM MANUAL 
"STOP" 


INTRODUCTION 

The FORTRAN IV version of the stellar opacity program, herewith described 
and referred to as the "STOP" program, is essentially a duplication of the IBM 
7090 FAP versions of the "IONIC" and "OPACITY" programs developed at the 
Los Alamos Scientific Laboratory by Dr. Arthur N. Cox. Other than the pro- 
gramming techniques employed, the only difference that exists in the present 
STOP version is that the "IONIC" and "OPACITY" codes have been combined 
into one program. In using the Los Alamos program to compute opacities below 
10 6 o K, the "IONIC" and "OPACITY" codes were run consecutively because the 
ionic method of computing occupation numbers was necessary for this low tem- 
perature region. After making a run with the ionic code to compute all equations 
of state data, the opacity code was entered to compute the absorption coefficients 
and opacities. For temperatures above 10 6o K, only a single rim was necessary 
with the "OPACITY" code, as this code employed another method of computing 
occupation numbers called the "Mayer Independent Electron" method which was 
applicable for high temperatures above 10 6o K. 

In the present STOP version, the user has an input option to determine which 
method to employ to compute occupation numbers. The two methods give the 
same values of opacity at temperatures near 10 6o K and this is normally the 
criterion for the user when making his runs. That is, for temperatures above 
10 6 o K, the straight opacity option would be used and for temperatures of 10 6 °K 
or lower, one would use the ionic option. 

The STOP program was designed primarily for use on the IBM 360/95 com- 
puter at the Goddard Institute for Space Studies, however, it may be used on any 
IBM 360 computer having a memory of 500 K bytes or larger, having a compiler 
which will accept the FORTRAN IV language . One should consult the USAGE 
section of this manual for further details on specifications. 

This manual has been written for the purpose of describing the methods used 
in the program and of showing the prospective user how to use the program. For 
those already acquainted with Cox’s works, accomplishing this is an easy task, 
however in attempting to write a manual which can be used by anyone wishing to 
compute stellar opacities, it was decided the following sections should be included 
in as clear and precise a manner as possible. 


a) General Description of STOP 


% 


This section gives a brief description of Cox's method of computing 
opacities and describes the general flow and usage of the program. It was 
designed for the person who is not acquainted with Cox's work and who wishes 
to find out what processes and techniques are employed. 

b) Users' Section 

This section describes in detail how to set up input to the program and how 
to interpret the results printed. The various options for intermediate output 
are described also. 

c) Cox's Continuous Opacity 

This section contains a detailed description of the methods used by Cox to 
compute opacity. Contained in this section is a description of the "Ionic" and 
"Mayer Independent Electron" methods for computing occupation numbers and 
other equations of state data. This includes the method of dominant potential 
and the equations used to compute the various scattering and absorption processes 
used in calculating the radiative opacity. A brief description of the conductive 
opacity and the procedure for adding this to the radiative opacity to give the con- 
tinuous opacity is also shown here. 

It should be noted that Cox's chapter in "Stellar Structures" was relied 
upon heavily in this manual and the notations used are the same as in the article. 

This project was begun at the suggestion of Dr. Albert Arking at the Goddard 
Institute for Space Studies, and it was under his guidance and direction that the 
present FORTRAN version of the program has been completed. Without his 
help in absorption coefficient and stellar opacity theory, understanding of the 
problem and ultimately completing the program would certainly not have come 
about. The invaluable help of John Stewart (1963) at the Los Alamos Scientific 
Laboratory cannot go without mention here also. His list of equations and flow 
charts of the Los Alamos FAP version of the program gave us much insite into 
many of the physical interpretations employed, and through his willingness to 
help and many phone conversations, we were able to complete this project. 
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CHAPTER I 

GENERAL DESCRIPTION OF STOP 


The present "STOP" program is a duplication of the IBM 7090 FAP, 
"OPACITY" and "IONIC" programs developed by Dr. Arthur N. Cox at the Los 
Alamos Scientific Laboratories. Its reason for existence is that it has been 
programmed entirely in FORTRAN. In describing the contents of "STOP," this 
section has been included to give the user a general idea of the processes used 
by Cox in his program and to also acquaint him with the general usage of "STOP," 
In a later section, a much more detailed description has included Cox's methods. 

Cox's method for computing stellar opacity employs the work of many men 
researching in the field of absorption coefficient theory. A list of references 
are given at the end of his chapter "Stellar Absorption Coefficients and Opacities" 
in Stellar Structures , however the basis for his method stems primarily from 
Harris Mayer in an unpublished Los Alamos Report, LA-647, (1947) entitled 
"Methods of Opacity Calculations." 

In this report, at the request of Edward Teller, Mayer studied the effect of 
line absorption on opacity for materials of extremely high temperature. He 
devised a method of computing the occupation of electrons in various energy 
states needed for calculating equation of state quantities for a one electron 
atom model, and calls this the "Mayer Independent Electron" method. He also 
describes another method of computing occupation numbers when temperatures 
are lower and a state of complete ionization does not necessarily exist and has 
referred to this as the "Ionic” method. 

Later, Cox, at the Los Alamos Scientific Laboratory studied the problem of 
calculating accurate monochromatic absorption coefficients and opacities over 
the whole range of temperatures and densities encountered in problems of stellar 
structure. He extended much of Mayer's work and developed a computer program 
which computes a Rosseland Mean Continuous Opacity given by 
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Employing both the "Mayer Independent Electron" method and the "Ionic” 
method to compute equation of state data, Cox's program, written in IBM 7090 
FAP language can be applied for stellar temperatures from 5 x 10 3o K to well 
over 10 7 o K. The "STOP" program can be used in this region as well. 

Stromgren (1932) transformed equation (1) in the form which has also been 
described by Chandrasekhar (1939) as 


where 
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X. S D. (U) 


Pj_ 
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kj (u) u 3 


N, ^ (u) u 3 
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and 


h v_ 

~ kT 


a. (u) is the cross section for Nj particles of matter density p. and mass 
fraction X. and is summed over all absorption and scattering processes. 

Mayer (1947) defined D in a different manner which is 


k i (u) Pj U 3 Nj or. (u) U 3 

D i ( u ) - A - A 


K ; (u) is the absorption coefficient per gram per nucleus for process, i , and 
A has inverse length units and is 


2 4 he 2 N 
3i3 mc kT 
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Nis the number of atoms per unit volume. The Stromgren D. (u)has the units of 
cm 2 /gm where the Mayer definition of D. (u) is dimensionless. 

The Mayer definition of D. (u) is the one used by Cox in his opacity calcula- 
tions and he has included the following processes in the computer program. 

a) Bound free absorption, D bf (u) 

b) Free Free absorption, D ff (u) 

c) Compton Electron Scattering, Ds(u) 

d) H~ bound free and free-free absorptions, DH bf (u)and DH f - f (u) 

e) H 2 + and H+H’ molecular absorption, D m ( u ). 

The induced emission factor is included in the definition of the weighting 
function w(u) and is given by 


15 u 7 e 2 


477 4 (e u - l) 


w(u) has been tabulated by Mayer (1947) however the program computes w(u) for 
a given u value. Below the plasma frequency cut off. Up, defined by 


Up 


477 1/2 M 1/2 

- (kT) \V) 


w(u) is defined to be zero and is limited in the program to fall in the limits 


Up < u < 30 . 


When u > 30, w(u) is so small and the D. (u)'s are so large that the value of 
the integral past this point can be neglected. The equations for the various 
D. (u) processes are described in Chapter III of this manual and are also de- 
scribed by Cox (1965). 
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Once the integral, I, is computed, then the radiative opacity can be computed* 
from 


K 


R 


_A_ 

Ip ' 


The units of k r using this definition is cm 2 /gm. 

In actual usage , the integration for I rarely needs to be computed beyond 
the maximum of the weighting functions at u ^7. An estimate of the remaining 
integral can be made by assuming the D. (u)'s will be no smaller at larger u 
values than at the edge where u = u e andD = D e , then 


AI 

n 




w(u) du 



e 


[S(«) - S(u e )] 


If this remainder is sufficiently small compared to the value of the integral at 
this point, it can be neglected. S(u) is the Stromgren Function and is tabulated 
as a function of u for use in the opacity code. 


r 

S(u) - l w(u) du 


In order to calculate several w(u) and D. (u) needed in the opacity integral, 
all the equation of state data including occupation numbers must be calculated 
from the users input of temperature, density, and composition mixture. Two 
methods are employed which were referred to earlier as the "Ionic” and "Mayer 
Independent Electron" methods. 

The "Ionic" method is the more accurate of the two and must be used for 
temperatures below 10 6o K, where ion formation is possible. For temperatures 
above 10 6o K, the "Mayer Independent Electron" method can be used with good 
accuracy and is the much faster of the two methods. Both of these methods are 
described in detail in Chapter in of this manual and by Mayer in (LA-647). 
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. Dr. Cox also described these two in great detail and has written some advantages 
and shortcomings in his article in "Stellar Structures" (1965). 

An important feature which Cox has included in the "Ionic" method gives the 
user of the program the ability to begin the calculation by specifying the tem- 
perature, composition, and either density, or electron pressure, or degeneracy 
of the gas. This is done using the method of "Dominant Potential" described 
by Chandrasekhar (1951). Since the method depends upon the degree of ionization 
of the mixture, it cannot be applied when very high temperatures are encountered. 
Therefore, when the "Mayer Independent Electron" method is employed, these 
options are not available and the only input is temperature, composition, and 
density. 

The calculation is much simpler in the temperature region 10 8o K up to above 
T = me 2 . Compton’s scattering by electrons is the most important process in 
determining the opacity in this region and Cox uses a method due to Sampson 
(1959). The opacity is given by 
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R 


N (etp) g C 


P 


G(T) 


where N (e+p) is the density of electrons and positrons, and G(T) is the Rosseland 
mean of the ratio of the effective Compton cross-section to o- Q , the Thomson cross 
section for scattering of radiation by electrons. 

At very high densities where electron degeneracy becomes important the 
conductive opacity must be considered. Cox uses the methods of Mestel (1950) 
for the conductive opacity. 


4ac / T 3 \ 

K c - 3 J • 

The conductive opacity is then added to the radiative opacity to give a total 
continuous opacity by 
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In using the program, firstly, a decision must be made as to which method 
of computing occupation numbers, etc. should be used. Since this can be based 
primarily upon the temperature, a general rule can be made. 

a. For temperatures below 10 6o K use the "Ionic" method. 

b. For temperatures above 10 6o K use the "Mayer Independent Electron" 
method. 

However, one should consult Cox's article (1965) for advantages and short- 
comings of this rule . 

For the "Mayer" method, one inputs: 

a. Temperature (°K, log l0 °K, kilovolts). 

b. Composition (number fractions, mass fractions). 

c. Density (gm/cm 3 ) 

and for the "Ionic" method, one has the following options: 

a. Temperature (°K, log! 0 °K, kilovolts), 

b. Composition (number fractions, mass fractions) and one of the following 

c. Density (gm/cm 3 ) 
or 

c. electron density (electrons/cm 3 ) 
or 

c. degeneracy parameter, v. 

Output options from the program give the user a wide range of flexibility. 
With no output options specified by the user, only the users input and final value 
of opacity is printed. These are given in table form if more than one case is run 
at one time. 

Various forms of intermediate output can be obtained through the output 
options and these include tables of energies, ionic abundances, pressures and the 
like. The "Users' Section" contains these options in detail. 
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CHAPTER II 


USERS SECTION 

The ‘'STOP' 1 program, which is the FORTRAN version of the Los Alamos 
Scientific Laboratories’ ’'OPACITY” and ” IONIC” programs, was written primarily 
for use on the IBM 360 computer although It is compatible with any FORTRAN 
system. It has been written completely In the FORTRAN IV language with the 
standard version In double precision. It has been written such that It Is 
compatible with any IBM 360 system having memory size of 500 K or more bytes. 

The ’’STOP" program has combined the "OPACITY” and "IONIC” codes Into one 
program which employs an Input option to determine which mode to run under. 
That is, the user supplies an input option specifying either an "OPACITY” 
run or an "IONIC ” run. The section. General Description of STOP , should be 
consulted for the calculations performed under each option. It suffices here 
to say that the "IONIC " option uses the "TuiIc” method for computing occupa- 
tion numbers and the "OPACITY” option uses the "Mayer independent electron 
method”. 

A complete package of the program consists of the source program 
(on punched cards) followed by the users input, and a data tape read by 
the program. 

The data tape consists of two files. Contained on the first file are the 
Moore element data tables and on the second file are the necessary tables for 
screening constants, gaunt factors, etc. Complete listings of these can be 
obtained from the authors . In running the program, the user would Include 
his data set behind the source program and place the data tape on logical 
unit 5. 

Included as Input are various output options which are explained below. 

If no output options are specified, standard output consisting of the users 
Input and tables of opacity, temperature, and density are given. 

A case is defined as one set of temperature, density, eta, or electron 
pressure, and composition mixture. The composition mixture can be further 
defined by the standard mixture ratios X, Y, Z If desired. Only X and Z 
are specified, however, as Y can be derived. 

As shown below, one data set can consist of several cases. That is, one 
temperature and several densities could be input and output would be a table 
of opacity versus temperature and density. Any combinations of temperature, 
density, electron pressure, or eta, and composition mixture can be run. 

Input consist of three categories: 

a) Mixture cards - contained on these cards are the output options. 

b) Temperature cards. 

c) Density, eta, or electron pressure cards. 
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INPUT 


'These are described as follows. 

All numbers should be right adjusted: 


columns 


CARD 1 


1 ~ 7 


19 - 20 

31 

32 

33 

34 

35 

k2 

50 


The format for this card is: 


variable 


MIXTURE 

NZ = no. of elements in mixture 

If non-zero, intermediate print-out for 
data tables 

If non-zero, intermediate print-out for 
opacity tab3.es 

If non-zero, print-out of initial conditions 

If non-zero, print-out of intermediate 
iterations for ionic code 

If non-zero, print-out of converged values 
and data passed to opacity code 

If zero: Mayer Independent electron method 

If non-zero: ionic method 

IFRAC = 1 - number fractions 

IFRAC = 2 - mass fractions 

IFRAC = 3 - basic set of mass fractions, X,Z 

IFRAC = If - basic set of number fractions, X,Z 

(A 6, 4X, 110, 10X, 1213., 18) 


Example : 

MIXTURE 13. 3.1111 3. 4 
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columns 


variable 


CARD 2 9 " 10 atomic number 

11 - 30 fractional amount as specified by IFRAC 

The format for this card is: (110, E20.8) 

Example: 

6 4.45E-4 

CARD 3 same as card 2 

CARD HZ 4- 1 same as card 2 

Mote : one card per element 

The following card is omitted if IFRAC = 1 or 2: 

CARD NZ + 2 1 - 16 X 

17-32 Z 


The format for this card is (2E16.8) 


Example: 

.5994 .0048 


CARD MZ + 3 1-5 

19 - 20 
29 - 30 


TEMP 

number of temperatures 
= 0 : temperature in kv 

= 1 : temperature as log 10 (temperature 

in degrees K) 

= 2 : temperature in degrees K 


The format for this card is (A 6, 4X, 2110 ) 


Example : 

TEMP 

1 

l 


CARD MZ + 4 

1 - 16 


TEMP (l) 


17 “ 52 


TEMP (2) 


33 - 48 


TEMP (3) 


49 - 64 


TEMP (4) 


65-80 


TEMP (5) 


The format for this card is ( 5 EI 6 . 8 ). Use as many cards as required, with 
at most five temperatures per card. 


lh 


Example: 

5-8 


CARD M2 + 5 


columns variable 

This card has to be one of the following three: 


option 1: 1-7 


DENSITY 


19 - 20 


number of densities 


30 

option 2: 1-3 


= 1 : density as gm/ cm 3 
= 2 : density as log 10 (density) 
gm/ cm 3 

ETA 


19 - 20 


number of etas 


30 


= 3 


option 3: 


1 - 10 


ELEC PRESS 


19 - 20 


number of el. pressures (j/CM 3 ) 


30 


= 4 


The format for this card is (A6, 4X, 2110) 


Example: 

DENSITY 2 2 


CARD N2 + 6 This card has the same format as card NZ + 4, depending on 

which option is selected. 

The above comprises a complete data set. Multiple runs may be submitted by 
repeating the above. At end of last set of data, there must be a card containing 
END in columns 1 - 5 . 


Example : 


MIXTURE 

2 


TEMP 

5-2 

DENSITY 

1.E2 


END 


1 11111 1 

1.000 

1 1 

1 1 


1 
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CHAPTER HI 


COX'S CONTINUOUS OPACITY 


In describing both of the occupation number methods as well as the opacity 
calculation as a whole, the following subscript notation has been used. The 
subscript, i , refers to the particular element present in the mixture and, j , 
refers to the number of electrons off element i. Therefore, if equation of state 
quantities are being discussed, say E^ , reference is to the total energy as- 
sociated with ion ij implying element i, with j electrons missing. The sub- 
script, k, refers to a particular state associated with either ion, i j , or 
element i. 

In the "Mayer Independent Electron" method, applicable for high temper- 
atures of 10 6o K or higher, (described by Mayer in LA 647, 1947, p. 26ff) all 
atoms are assumed to be completely ionized and the electron wave functions 
approach those of a free electron and are therefore independent of the positions 
of the nucleus. In this case, only the position of an electron in the field of the 
nearest nucleus needs to be computed. 

Therefore, the electron state energy in this method need only to be re- 
ferred to by E. k . That is, in the "Mayer Independent Electron" method one 
computes equation of state quantities using only the subscripts i and k, which 
refer to a particular nucleus i, and state k, associated with this nucleus. 

In the "Ionic" method, used primarily for temperatures below 10 6o K, the 
assumption is that an electron is in a state of sufficiently low energy that it will 
be temporarily bound to a given nucleus. For each nucleus which may have 
several electrons bound at the same time, the one electron wave functions will 
depend upon all the bound electrons for each nucleus. Therefore, one must speak 
of the number of ions of each type. This has the advantage of taking into ac- 
count the interactions of the free electrons and ions. 

Therefore in the "Ionic" method all three subscripts i, j, and k are used 
to denote, say the electron state energy E. jk , which refers to the energy of an 
electron in state k of ion i j . 

For example, in referring to Cox's equation (143) for bound-free absorption, 
in his article in Stellar Structures, the sum goes over all states of all the nuclei 
in the mixture with ionization energy less than (u • kT) and equation (143) would 
take the following forms for each method: 
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MAYER INDEPENDENT ELECTRON METHOD 

For a given composition, temperature and density, all electrons are as- 
sumed to be ionized from the atoms. The number of free electrons per average 
atom is then 



i 


Using the Fermi-Dirac distribution, the number of electrons in bound state k of 
element i is 


X 


ik 


k 


'ik 


kT 


1 


where 

g ik is the degeneracy of state k and 

E ik is the interaction energy of electron in state k with the nucleus. 
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For this first approximation, the total number of bound electrons per 
average atom in the system is 


XB = 



X. X 
1 


ik 


and a new number of free electrons is 


y ~ 



These equations are repeated until successive y values differ by a predetermined 
amount. 

The energy, E. k> of the bound electrons, which includes perturbations by 
other bound electrons and continuum depression due to free electrons is due to 
Mayer (1947) and is given by 


E 


ik 


Ryd 


Z *k 2 


Vi_ 
r . 


3 - 


r . 2 Z * 2 
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+ 5 
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X. y. 
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v r i 


with the effective charge given by 



where the quantity 


S 


k 


z 


*2 
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is tabulated by Mayer (1947) r k is the average radius of state k in Bohr radii. 
The quantity 



is limit to be £ 1. If > 1, it should be set = 1, for calculation. 
The quantity on the right in the energy equation, 


.3 

5 



1 


X. 


y 


2 

i 


y r i 


is the negative of the potential energy of the free electrons averaged over all 
ion spheres and is usually called E 0 . The screening constants cr kk , are tabulated 
by Mayer also, r the average radius of a sphere which contains enough elec- 
tronic charge to neutralize the ionic charge of atom i, is 


r . 



(Bohr Radii) 


where V, the volume per average atom is 


V = 


M c 

PN 0 



P v o 


and a 0 is the Bohr Radius. 

This method for obtaining occupation numbers converges rapidly except when 
atoms are nearly neutral or for temperatures greater than approximately 10 6o K. 
However, it does not work for stellar atmosphere temperatures and densities 
because the energies needed in the occupation number expression is very sensi- 
tive to the presence of other bound electrons in the same atom. 
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The degeneracy parameter, 97 , for free electrons is obtained in the Los 
Alamos code by an iteration procedure using a table of 7 / vs. Fermi-Dirac integrals 
Fj 2 ( 17 ) depending upon the state of the degeneracy. This procedure is outlined 
below. 

For rj <-28, no iteration is necessary and 


V 


ln [ F 1/2 (v 



2 



For 97 > -28 


A first guess is made of 


277 2 

^ 1/2 ( 7 ?) (kT ) 3,/2 V ^ 


after which a table of 97 vs. F 1/2 (97) is entered to interpolate (quadratic) to ob- 
tain 77 .. Using this value of 97, a new value of F 1/2 ( 77 ) is computed from one of 
the following: 



= (O.OO 877 3 + 0. I 8797 2 + 0. 525t7 + 0.678) 

+ (-0.012t ) 2 + 0.01977 - 0.009) 1 77 1 ; 

- 0. 8062c 77 - 0.3133e 27? + 0.1569e 377 ; 


rj > 30 
4 < rj < 30 

-2 < 77 < 4 

-28 < r; < -2 


28 


and the table is re-entered to obtain a new 77. + This process is repeated until 
a consistent 77, or 


Vi - Vi-i 


Vi 


< e„ 


is obtained. The program has set e v = 10 and for the first iteration 77. - 2 

x 10 4 . Also an average value for the number of bound electrons in each energy- 
state for each iteration given by 
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i k 


i k 


+ x. 


1 k . . 
old 


At this point when a consistent number of free electrons, y, occupation 
numbers, X. k , etc. exists the equation of state quantities can be calculated. 

The matter energy equation of state is computed by adding the energy of 
ionization and excitation to the kinetic energy of the free electrons. If the 
densities are high, electrostatic corrections for interactions between all the 
free charged particles must be made to both the energies and pressures. The 
energy expression given below by Cox differs from that given previously be- 
cause the interaction energy between free electrons and the nucleons as well 
as the potential energy of the free electron cloud are included. 

This energy for nucleus i is 


E. 

1 Ryd 



7 **2 
ik 




9 yi 2 2 

"5 r. + 2 ( kT )Ryd. 


where n k is the principal quantum number of state k. Z. k is the effective 
charge seen by the electron in state k and is 


7 ** 
ik 





Cr 


kk 
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The sum is over all states which have energy less than that in state k. 
cr kk are the Slater screening constants and as before, are in tabular form for 
use in the code. 

The total matter energy is given by Mayer and is 


"Ryd 


L 


X. E. 


V(kT), s/2 


f R y d 


2v‘ 


3/2 


(V) 


where F 3/2 (rj) is the Fermi-Dirac integral of order 3/2. The total pressure 
equation which is from Mayer also is 


F (Ryd/B.Ryd 3 ) 



i 


with the electron pressure given by 


P 


(kT) Kyd 

V 


2 (MV/d 

Y — 

3 2n 2 


A i 

F 3/2 (V) - 5 (kT) Ryd 






At this point, all quantities needed to complete the opacity integral 


I = 



w(u) 



Dj (u) 


du 
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i 


, are tabulated and pertinent quantities are sorted on increasingly negative energy 
E. k , and u is computed from 


u 


E. 


ik 


(kT) 


R y d 


Since the integrand for I is discontinuous at each u value or edge, integra- 
tion must be performed between each u value separately and then added together 
.to form the total integral I. Integration between each edge is done using 
Simpson's rule. 


IONIC METHOD 

For temperatures in the range of 5 x 10 3o K to 10 6o K, Cox has developed the 
following method, referred to as the "Ionic" method, for calculating occupation 
numbers. 

One starts this method by calculating the equation of state of a given mixture, 
temperature, T (°K) and degeneracy parameter, 17. The form of the Saha Equation 
applicable to any degree of electron degeneracy is 


i.i + l 
i . j 


B. 


1 ■ j 


+ 1 


B. . 
1 . j 


exp 


( J + 1) e 2 

E 0 

kT kT 


where 

x. - no. fraction for ion i • j , i = element, j = no. of electrons 

missing 

B - partition function of ion i, j 

Xi j - energy needed to ionize ion i, j to i, j + 1 

( j + 1) e 2 /R D - correction to ionization energy for interaction between 
charged particles 
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V 


- degeneracy parameter corrected by E 0 /kT 


E 0 - potential energy of the free electrons-averaged over all ion 

spheres. Enables one to ignore interaction energy of the 
free electron with themselves and with other ions 

R d - Debye Radius. 

The Debye Radius considering all charged particles is 



The partition function for ion i j, summed over all excited states from the 
ground state to principal quant. no.n raax = (j + 1 )r B Rad 



gj . k statistical weight, if not tabulated, then either 2n 2 or 2(21 + 1) 
e . . k excitation energy alone the ground state 
y i electrons off element i 

r. average radius for element = (3/47T y. /y v) 1/3 /a 0 
V volume /average element 

a 0 Bohr radius 
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Ryd hydrogen atom ionization energy - 136 ev. 

5 = r. 2 Z* 2 Bohr rad. 2 tabulated by Mayer ion sphere model. 

n k 

e . . k is either tabulated by Moore energy level data or computed using 
screening constants. 


Z 


* 

i jk 


Z. 



v k ^kk 


and 


i j k 


7 * 2 

ijk 

n ijk 


XiJ 


e . k may include states where more than one electron is excited above 
the ground level and are tabulated from Moore "Atomic Energy Level Data," 

NBS Circ. 467, 1, 2, 1949 and 1952. 

If level data does not exist for a particular state, screening constants 
given by Mayer and Karzas are used with the Bohr model of the atom to compute 
these data. Hydrogen like statistical weights are used, that is gj . k = 2n 2 , 

(n >6), or 2(21 + 1), (n <6). 

In cases of low temperatures where most atoms are neutral, the perturba- 
tion of the energy levels is due principally to neutral atoms. Therefore, a 
partition function is cut off when the orbit of the excited level has a radius 
greater than the mean separation between nuclei. 

The procedure for solving the ionization equation is as follows. Once 
initial values for the ionization equation is made, one solves the two polynomials 


no. of nuclei of ele ment i 

P . - ~ y _ 

1 n 


X. 
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X XX 

no. of electrons from element i ii _ i2 A ii 

= jt: = X7 +2 


Xu x i0 


+ z 


l il 


i z- 1 


Therefore, y., the number of electrons from element i, is 


yi 



free electrons from element i 
no. of nuclei of element i 


and the total no. of free electrons per average element is 


y 



L y . 

i J i 


The electron density is then computed from 
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47r(2nkT) 3/2 
- 3 - 


F 1/2(v) 


and the average atom density 
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N 0 
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P-c 


where, p, the matter density is 


P ~ 
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PRECEDING FACE BLANK NOT FILMED. 

If the case is non-degenerate 
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The occupation numbers can then be calculated from 
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Corresponding effective charges and energies for each level along with ionic 
abundances are tabulated for future use in the absorption coefficient calculation. 


METHOD OF DOMINANT POTENTIAL 

Cox has developed a method using the method of dominant potential de- 
scribed by Chandrasekhar (1951) whereby one can input either temperature and 
density or temperature and the state of degeneracy, 77 , or electron pressure, 

Pe for the computer program. 


The method, as described by Chandrasekhar is as follows: 


Letting x (n) » N e , N (n) and N n 1 be the ionization energy, the electron 
density, and atom density with (n) and (n - 1 ) electrons, then 


N ( n ~ 1 } 
]\j(n) 


; (n_1) (2n nkTl 3/2 


r(n) 


- - ( n ) A.n 

e A / 
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where g ( ” 1 ) and g (n) are statistical weights. Then the state of ionization can 
be inferred from the value of a dominant potential 0 defined by the equation 


1/,/kT - JL ( JE nkT) 3/ f 

e _ Nc h 3 


because, if x (n) = 0 then apart from the statistical weights, the number of atoms 
with n electrons will equal the number of atoms with (n - 1) electrons. In a 
first approximation one can say that an atom with nuclear charges, z, will be 
mostly (z - n) fold ionized if y (n) - 0 . 

The method, as employed by Cox, is described below, and after one of the 
following methods are completed, the ionization equation is computed as be- 
fore assuming that the temperature 6 - kT, degeneracy rj, and the density p are 
all self consistent. 

Temperature, Density Input 


First the atom density, N t , is computed from 


N 


* (a tom/'cm^) 


No/; _ Ne 

~ y 


where No is Avogadros number, p c the cold mean molecular weight = I X. A. 
and y the total number of free electrons per average atom. 

The dominant potential equation 


^/kT = _L ( 277 n kT )- 


3/2 


Ne 


is written as follows to calculate the known quantities. 


0 

lny + kT 


In 


(d 3/2 \ ( 2( 2tt n) 3/2 \ 

V N t / + ln \ h 3 / 


DPC 
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The unkno\vn quantity on the left, ( In y + 0/ kT), will be referred to as the dom- 
inant potential constant DPC. This equation is iterated using the known ioniza- 
tion energies for each ion, i j, in the mixture to get a consistent pair of v and ,/; 
for the given temperature and atom density. 

The procedure is as follows: 

Using a starting guess of 0 = 5 • (kT) kv, compute for all elements in the 
mixture, a number of electrons off element i , y . from the following 


Vi - 

where is the number of electrons 
than the ionization potential, x'. ■ to 
pute the step y 1 from 


j + e kT 


1 


k v 


off element i for the ion when ,■ is just less 
ionize the atom to ion i, j • 1. Then com- 


y 


i 



X 


i Vi 


and a new value for 0 from DPC. That is 


0 


n e w 


(kT) kv 


DPC - 


In 



where y y'. Then compute a new' set of yds, y and and keep iterating 
until succeeding y 1 values converge to some predetermined amount . The 
code has =0.1. 

Once y is known, then the electron density, Ne, and electron pressure, Pc, 
can be calculated using Fermi-Dirac statistics as follows. First compute the 
electron density, N c 


Ne N t • y 
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and the Fermi-Dirac integral of order 1/2 in the non-relativistic case 


F l/2 (V) 


h 3 Ne 

4tt( 2n kT) 3/2 


Knowing this, existing tables of 77, F 1/2 (77) , F 3/2 (77) can be entered to com- 
pute the other two and compute the non-relativistic electron pressure from 


Pe 


Ne(kT) 2/3 


F 3/2 ( V ) 

F 1/2 C 7 ?) 


477 ( 2n) 3/2 
h 3 


(kT) s/2 


3 F 3/2 



At this point, the ionization equation can be solved assuming a consistent set of 
Ne, y, Pe, 77 F 1/2 (77) and F 3/2 (77) for the given input of temperature, (kT), and 
density, p. 

(Temperature, degeneracy) or (Temperature, electron pressure) Input 

This method is started by first computing either the electron pressure, Pe, 
or the Fermi-Dirac integrals depending upon the input for the non-relativistic 
case from 


Pe 


2_ 4tt( 2n) 3/2 

3 F 3/2 C 7 ?) .3 


(kT) s/2 


F 3/2 (77) and F 1/2 (77) can be obtained in a table from the inputted, 77. If Pe was 
input the above equation is solved for F 3/2 (77) obtaining 77 and F 1/2 (77) from 
tables vs. 2/3 F 3/2 (77). 

Then with Pe, 77, Fj /2 (77) and 2/3 F 3/2 (77) known, the electron density is 
computed from 


Pe F i/2( 7 ?) 

(kT) ^ F 3/2 (v) 
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Ne 


and the dominant potential <p from 


e4 /kT = 2 

Ne h 3 


The number of electrons off element i, can be computed from 


Vi 


y + e 


'/'“Xjj/ (kT) 


kv 


where j is chosen as before. 

Now the total number of free electrons from an average atom in the mixture, 
y , the atom density Nt , and the matter density p can be computed from the 
previous equations 


y - 



Nt 


Ne 

y 


Nt Mc 

p ~ No 


Now, as in the temperature, density input, one has a consistent set of all 
quantities needed to begin the ionization equation calculation. 

After the ionization equation is solved, the same procedure is followed as 
in the "Mayer Independent Electron" method to solve the opacity integral, I. 


to 


First all quantities are tabulated and sorted on increasingly negative E . 
u is computed from 1 1 k 


u 



(kT) 


kv 


and the integration is performed between each edge using the trapezoidal rule. 


THE PROCESSES 

Bound -free Absorption, D bf 

Cox writes the following equation for D bf (u) (equation 143) 


(a) 
b 



T X. X. . X. 

/ ■ 1 1 J 1 J k 


477 4 n 2 e 8 
h 4 n k 5 (kT) 5 


f ^b f 


where K bf is the absorption coeff/gm for ion i j and is summed over all energy 
states, K per ion, that is 


K 


ff 


L X. X . . X. .. 

1 i J i 1 k 


bf . .. 
i Jk 


N_ 

P 


The cross section for bound free absorptions o- given by Cox (equation 23) is 


2*n' 

3/3 ch 6 n s v 3 qb < E « 


X. is the number fraction of element i, X. is the abundance of ion i j 
(element i, ionized j times) such that X. X. is the number fraction of ion ij, 


4i 


and X. . k is the occupation no. or probability that an electron exists in state k 
of ion i j . 

5{C ^ 

Z k is the effective nuclear charge seen by an electron in state k consider- 
ing screening to both bound and free electrons, n is the principal quantum num- 
ber, v the frequency (sec -1 ) and the other atomic constants have their usual 
meaning. 

q b f is the correction to the cross section for the availability of the final state, 
and is 


q b f 


■ e n -u-(E ijk /k T ) + ^ 


t) defines the matter state of degeneracy and will be referred to as the 
degeneracy parameter from here on, 


u ' kT 


and according to Mayer's A definition 


2 4 he 2 N 

3/3 me kT 


q u is the gaunt factor making cr bf the non-classical cross section. In comparing 
Cox's published equation with that in the code, or 


D bf 


yh z *' 4 yxn' 

( kT > R 2 yd n 1 5< qbf i ,gbf i' 
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a discrepancy of the factor 


47 r 4 m 2 e 8 


is discovered, however this is exactly the energy of one Rydberg squared and 
cancels when using the temperature in Rydbergs. 

The "STOP" uses a set of tables of computed bound-free gaunt factors ob- 
tained from John Stewart at the Los Alamos Scientific Laboratory which uses 
a polynomial fit approximation developed there of the Karzas and Latter (1958a, 
b, 1961) bound-free gaunt factors. 

Cox states, "The bound-free gaunt factors depend on the initial state quan- 
tum numbers and photon energy in excess of the ionization energy. All extensive 
opacity calculations assume that bound-free transitions occur in hydrogen-like 
atoms where the electric fields are Coulomb, and in this case, the free electron 
energy parameter on which g bf depends can be scaled by the effective nuclear 
charge Z ." 

The most recent and most accurate hydrogen-like gaunt factors are those 
by Karzas and Latter (1958a, 1961, 1958b). These are functions of 


Y 2 

i 


(u-u e ) 


where 

u e = u at the present edge. 

The effective charge > is defined by 


Z 


* ' 


- ( E , • ' E„) 


1/2 
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where E ( < represents the electron energy in state 1 ' and E 0 is the potential 
energy of the free electrons averaged over all ion spheres. This shift in the 
energy zero point enables one to ignore the interaction energy of the free elec- 
trons with themselves and with the ions which are assumed to have bound elec- 
trons with orbits small compared to 7 S *7. is the average radius of a sphere 
which contains enough electronic charge to neutralize the ionic charge of atom 
i . E 0 is defined by 


E 


o 



y ; is the number of free electrons off element i, (and is taken as an ap- 
proximate effective nuclear charge seen in state k in later calculation) 




and 


yi 3 

7^ V 


) 7 - 


whose V is the volume per average atom 


Me 


z>- 


PNq PNq 


a 0 is the Bohr radius, N 0 is Avogadros number and p c is the mean molecular 
weight, p, is the density in gm/cm 3 . 


bb 


The summation is over all edges which have ionization energy less than 

(kT). 


Free -Free Absorption D ff 


Cox has derived the following definition for D f f (u) which is used in the code. 


277 2 me 4 y i 
i h 2 kT ln 


with K ff the absorption coeff/gm and may be summed overall particles in the 
mixture. 


( 1 + ^) g ff 


D ff (u) 


k ff 


L 


K 


f f 


N 

> X. X. . O-., . - 
/ . 1 1 J f f ij P 


where the sum goes over all ionization stages, j, for each element i, N is the 
particle density and p the matter density: The cross section cr. is given by Cox 
as 


CT . . 

1 J 


2 s 7T 2 Z* 2 e 6 kT 

3y3h 4 cz7 3 ln 


1 + e 77 


gff 


The free-free gaunt factors g ff are those by Karzas and Latter (1958a, b, 
1961) and have been averaged over the Fermi-Dirac distribution and corrected 
by Green (1958, 1960) for screening of the nuclear charge and availability of a 
final state. The Karzas and Latter data is in tabular form and the Green correc- 
tions are applied in the program. 


Green writes. 


gp D ( u ’ ^ 2 - ^ = i 1 



- 4 


log 


(l + e v ) (l + e Tr ' 2 ~ 2 ' y A) 

(l + e 77-11 ) (l + e a ~ 2y Z 1 ) + gFD^T - > u > r r 


l > ' / 
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. q 

where g F ° u , the un-normalized gaunt factor including screening, has been writ- 
ten in terms of the tabulated Karzas and Latter gaunt factors g FD and g 0 („;) , 
the asymptotic value of the unscreened hydrogenic gaunt factor for E = 0. 
g 0 (a).) has been tabulated by Karzas and Latter also however the program uses 
a polynomial approximation to compute g 0 as a function of oj.. is de- 
fined by 


u 



where u has its usual meaning and 

Z 2 

y? = rn 


t ) is the degeneracy parameter and 1 is a dimensionless length such that the 
screened potential is given by 


V 


2e~ , °/ 1 

P 


Results obtained from the polynomial for g 0 O) we checked against Green's 
tabular values and agreement was only to one or two figures. Green states "In 
the classical case, 1 is just the Debye length divided by a 0 /Z. At higher den- 
sities, insofar as the one-electron approximation can be applied, the potential 
that the electron sees must go into the Fermi-Thomas potential. Kidder 
(UCRL 5267-T, June 1958) discusses the behavior of this potential and points 
out the circumstances under which it may be approximated by the above potential, 
V, with a value for la 0 /Z somewhat different from the Debye length." 

The Los Alamos code uses the following for 1, 


r 1 


R D /a 0 




Finally, the screened gaunt factor used in the code is calculated under 
conditions of electron degeneracy, 77, as 


iff [y 2 > u - v) 


fr i F 5 D (' y 2 ’ u - v) 

21 1/2 C 7 ?) 


where I 1/2 (7?) is the Fermi Dirac integral 


1 1/2 C 7 ?) 


In actual usage, three equations have 
various degrees of degeneracy. 

Cox states "for free-free absorption the nuclear charge seen by the free 
electron is taken as yi , and even if data for individual ions are known, the sum 
is only one term per element. This is possible because free-free absorption 
is never important when there is only partial ionization." 

The factor, 2tt 2 me 4 /h 2 , here as before, cancels when using (kT) in Rydbergs 
and the form of the equation in the code is 


f 

«*0 


& 

1 + e x_7? 


dx . 


been used in the code depending upon 


D 


f f 



In (l + e 71 ) g ff (u, y . 2 /kT, 77) 


Scattering Ds(u) 


Cox writes the following equation for D s (u) 


D s <U) 


-° u3 


A 



1 


u 3 V3 

777 ^^)^ 



h 2 

2tt 2 me 4 


G(u) . 
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G(u) is the ratio of effective electron-scattering cross-section, g f f , to the 
Thomson cross section for scattering of radiation by electrons, a 0 , given by 
Sampson (1959) that is 


G(u) 


cr 


ef f 


cr 


0 


where 



Cox writes 


Cr 

e 



V -N 

X i Vi CT eff p 



1 


x i y i a o G ( u > 7T 


where X t are the number fractions of element i, y ; the number of free elec- 
trons from element i, N the particle density and pthe matter density, u and A 
have these usual forms defined in equations 65 and 67 in his chapter in "Stellar 
Structures." 

cr g is the absorption coefficient for electron scattering/gm and is summed 
over all elements in the mixture. The induced emission correction (l - e“ u ) for 
pure absorption processes should not be applied to free-electron scattering and 
is therefore put into Ds(u) to cancel out the factor in the weighting function w(u). 

G(u), which is the same as G (t , t ') in equation (31) in Sampson (1959) can 
be represented by the following expression (as derived by Cox) 


G (u) - 1 + B(T) ( c ° + c 2 u + C. 


u 2 + c, 


3 X 

u + 


u 4 + c, 


u S ) 
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where 


B(T) = 1 + 1.875T + 0.8203125T 2 - 0.3076617188T 3 + 0.367270229T 4 - 0.51549911T 5 
c 0 = +(2T + 8.75T 2 + 14.76525T 3 + 6.7675781T 4 - 2.7493286T 5 ) 
e j = -(3.2T | 26.6T 2 + 122.85T 3 + 351.91406T 4 + 549.86572T 5 ) 
c 2 = +(10.5T 2 + 141.4875T 3 + 1011.8883T 4 + 4600.3966T 5 ) 
c 3 = -(31.4714285T 3 + 620 .3160 7T 4 + 6341.6191T 5 ) 
c 4 = +(88.114286T 4 + 2380.7571T 5 ) 
c 5 = -(234.85714T 5 ). 

T is in units of me 2 and is T = (kT) k v /510.974 kv/mc 2 . This is the form used 
in the code. 

The equation used in the program for D s is 

u 3 f2ne 2 \ 3 

D s< u ) = (kT) G(u) 


where the factor 


h 2 

277 2 me 4 


is, of course, included in the temperature conversion for (kT)in Rydbergs. 
The total number of free electrons per average atom, y, is 


y 



X. y. 

1 J l 


1 


^9 


H“ Bound- Free Absorption D H ( - J 


Cox writes in equation (147) 


D h- u) 

H bf 


K„- pu 3 

H bf 

A~ 


V 

a «i, v x » x »- 


3^3 me 


(kT) u 3 


where, using the Saha ionization equation at the temperature and density of in- 
terest with the known ionization energy of .75 ev for the H" ion gives the fol- 
lowing for the ratio of the H" ion abundance per neutral hydrogen atom, 


X H - 

YT 11 4. 158 X 10” 10 P tfS/2 1.7260 

A„ O e 


where r = 5040/T°K and other units are in cgs. 

The bound-free H - cross sections, <r H - , used are published by Gelfman 
(1962), and here as with the H7 f absorption, the H" ion is not allowed to exist 
above T = 2V = 23, 210°K. 

The tables, have Avogadros Number, N 0 , included in them, or N 0 ct h _ and 
are tabulated as a function 


(( kT > R y d • U " E o) 


The form of the equation used in the code is as follows. 



N 0 °-h 


b f 


V 3 

v u x hV 


£_ 
V- c 


a 0 V 
AV 


Equivalence between this and Cox's published equation can be seen from the 
following. 
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The average atom density, N, can be written as 


N 


Ne 

y 


P 


N 


o 


Pc 


where Ne is the electron density, y the average number of free electrons per 
average atom, p, the density, and pc , the cold gram molecular weight. The 
quantity a 0 V/AV is just A" 1 in centimeters, and therefore 


N 


o 


Pc 



3V3 me 
2 4 he 2 


(kT) 


H~ Free-Free Absorption P H ( - U ^ 

Cox writes in equation 148 in "Stellar Structures" 


D 


f f 


Kh- pu 3 

H ff 

A 



P. X. 


X H u ; 


3/3 me 
2 4 he 2 


(kT) 


where the absorption coefficient per gram, K H - , is tabulated by Ohmura and 
Ohmara (1960, 1961). The quantity tabulated in f their article is the continuous 
absorption coefficient of the negative hydrogen ion per neutral hydrogen atom 
and per unit electron pressure for free-free transitions after allowing for 
stimulated emission factor. 

Cox has tabulated the following quantity for use in the code 


(Ak 2 ) 3 K„- ( N 0 


p , V 


where (Ak 2 ) = u(kT) RyJ and is tabulated as a function of (kT) kv and (Ak 2 ). 
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The form of the equation in the code is then 



The equivalence of this equation with that by Cox is of course seen as before 
by knowing that 


3/3 me 
2 4 a 2 

e 


(kT) 



H 2 + Absorption D H + (u) and H+H Abs D H+ (u)H 
Cox writes (equation 149) for molecular effects 


D (u) 




DH 2 + 


+ D 


H + H 


a) DH 2 + 


dh 2 + 


V 


3^3 me (kT) 
2 4 he 2 N t 


A h 2 + is tabulated by Bates (1952) and are only applied in the temperature range 
(2500 - 12,000°K). The code has 


d h 


+ 

e 


K^+u 3 X h 2 XB h X(H, 
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Since 


/> 



and 



The tabulated values of K H +are scaled by 10 39 . 
b) ^h + h 

Here, Cox writes 


D 


H + H 




X H X H 0 ^ 


3V3 mc(kT) 
2 4 he 2 N t 


where A H+H is tabulated by Zwaan (1962) and is applied in the temperature 
range (3000 - 8000°K). The code has 


D 


H + H 


10 logA H + H U 3 


( p - p .) 


2 Y 2 
A H 


a 0 V 
AV 


xb h 2 


TEMPERATURES ABOVE 10 8o K 

In the temperature region 10 8o K up to above kT = me 2 , Compton scattering 
by electrons is the most important process in determining the opacity (Sampson, 
1959), and Cox uses the following for the opacity resulting from Compton scat- 
tering alone. 


K 


c 


N (e+p) CT 0 


P 


G(T) 


( 1 ) 


This equation neglects absorption and N ( e + p > = density of electrons and positrons 
G(T) is 


G(T) 
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and is the Rosseland mean of G(u, T' ). G(u, T' )is the ratio of the effective 
Compton cross-section, defined by Sampson, to o - 0) the Thomson cross-section 
for scattering of radiation by electrons. T' is in units of me 2 and u = hzVkT 


a 


o 



6.65209 x 10" 2S cm 2 


For the temperature range 20 kv < kT £ 125 kv, Sampson derived the follow- 
ing polynomial fit for G(T), T is in units of kv 


G(T) = -0.13887 + 4. 9871(T)~ 1/2 - 5.9479(T)" 1 - 2.362(T)" 3/2 

Equation (1) can be written as 


K 


C 


N( e+P) a 0 N 0 _ 

Nq — G(T) 


0.40077 


N 


( e + p ) 


PN 0 


G(T) 


From equation 132 in Cox’s chapter 


N 


e 


PN 0 



where he has neglected any positron concentration and therefore obtains the 
following equation for k c 


K 


C 



( 2 ) 


where 

G(T) = -0.055655 + 1.9987T" l/2 - 2.3837T" 1 - 0.94662T" 372 


Cox uses equation (2) down to kT = 14 kv which gives rise to a small difference 
in the computation of G(T) and therefore to the radiative opacity, k c , as well. 
G(T) was checked at kT = 14 kv to compare the result using the above poly- 
nomial with the tabulated value published by Sampson. Sampson's value was 
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* 

obtained by numerically integrating the integral representation for G(T). 


^^^Sampson 

~ 0.7279 

(kT = 

14 kv) 

G(T) Cox = 

0.7241 

(kT = 

14 kv) 


ELECTRON CONDUCTIVE OPACITY 

At very high densities where electron degeneracy becomes important, one 
must consider electron conductions with radiation flow for transportation of 
energy. 

The conductive flux can be expressed in the form following Mestel, (1950) 


T 3 dT _ dT 
K c p dx ~ v c dx 

(Cox, p. 222) from one of the three cases: 

1) rj < -4 (non-degenerate gas) 


_ 4ac 

F c = - T~ 

where the conductivity v c is obtained 


128 mk 5 T 4 e 7 ? 


h 3 e 4 pN 


°L 


y . X . 0 . 

J i i i 


Mr 


where 


1 2 

@i - 2 1 - cos co. 
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and oj i is given by 


fez 173 ] 

CO. ~ -rji - 

1 Vi 173 


and for this non-degenerate case 


ez 


1/3 - 


0.58926 e v/ 3 


2 ) 17 > 8 (high degeneracies) 


167T 2 mk s T 4 (t ? 3/2 + -Xt ] - 172 ) (l + 9.376 (1 


-4/3 


T ) 372 + -g" T) 


~ 1/21 


9h 3 e 4 pN 


-L 


y . 2 X. 0. 

J l 11 




(1 + X 2 ) 


and 


f 2 - 


1.02715 x 10 


■6 


2/3 


877 m 3 c 3 
3 h 3 N„ 


10 e 


1.02715 


where 0 j and d i are given as before with the exception, 

i) [$Z 1/3 ] is obtained by interpolation from a table by (Cox, p. 223) of con- 
ductivity integrals for the interval -4 5 77 < 30 


ii) [<9Z 1 7 3 ] = .84 753, 77 >30 
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3) (-4<t]< 8) (partial degeneracy) 
Cox has tabulated the quantity 


V^i 2 x i e i 

v ‘ p 


for this case. Therefore one calculates©., as specified in 2) and divides out 
the known quantities to obtain v c . 

Now k e is obtained from the flux equation 

4ac / T 3 \ 

"e = 3 [pvj 

a = radiation constant = 5.6692 xlO’ 5 erg cm -3 deg’ 4 
c = velocity of light = 2.997929 x 10 10 cm sec’ 1 

The conductive opacity is then added in the following manner to the radiative 
opacity to give a total continuous opacity * T by 


1_ = J_ J_ 

C T K R K c 
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